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7 -ray spectra for positron annihilation with the core and valence electrons of the noble gas atoms 
Ar, Kr and Xe is calculated within the framework of diagrammatic many-body theory. The effect of 
positron-atom and short-range positron-electron correlations on the annihilation process is examined 
in detail. Short-range correlations, which are described through non-local corrections to the vertex 
of the annihilation amplitude, are found to significantly enhance the spectra for annihilation on the 
core orbitals. For Ar, Kr and Xe, the core contributions to the annihilation rate are found to be 
0.55%, 1.5% and 2.2% respectively, their small values reflecting the difficulty for the positron to probe 
distances close to the nucleus. Importantly however, the core subshells have a broad momentum 
distribution and markedly contribute to the annihilation spectra at Doppler energy shifts > 3 keV, 
and even dominate the spectra of Kr and Xe at shifts > 5 keV. Their inclusion brings the theoretical 
spectra into excellent agreement with the experimental 7 -spectra across the full range of Doppler 
energy shifts. Additionally, the theory enables the calculation of the ‘exact’ vertex enhancement 
factors 7 'ni for individual core and valence subshells n£. They are found to follow a simple and 
physically motivated scaling with the subshell ionization energy I„e: 7 n^ = 1 + sjA/I n i + (B/Ke )^, 
where A, B and /3 are positive constants. It is demonstrated that these factors can be incorporated 
in simple independent-particle-model calculations to successfully reconstruct the true many-body 
annihilation 7 -spectra. This formula can be used to determine the enhancement factors for positron 
annihilation with core electrons of atoms across the periodic table, and with the localized atomic-like 
core electrons of condensed matter systems. 

PACS numbers: 78.70.Bj, 34.80.Pa, 34.80.-i, 34.8.Uv 


I. INTRODUCTION 

Low-energy positrons annihilate predominantly on the 
outermost valence electrons in atoms, molecules and in 
condensed matter systems [1, 2]. Small fractions of 
positrons can, however, tunnel through the repulsive nu¬ 
clear potential and annihilate on the core electrons [3]. 
The 7-rays produced in the annihilation events give rise 
to a Doppler-broadened energy spectrum that is charac¬ 
teristic of the electronic states involved, e.g., annihilation 
on the tightly-bound core electrons results in a distinct 
signature in the high-energy region of the spectrum. A 
dramatic role is played in the annihilation process by 
short-range electron-positron and positron-atom correla¬ 
tions: they significantly enhance annihilation rates [4, 5] 
and alter the shape and magnitude of the annihilation 
7-spectra [6-10]. 

In this work we examine, systematically and in de¬ 
tail, the effect of these correlations on the 7-spectra for 
positrons annihilating on individual core and valence sub¬ 
shells in many-electron atoms. Specifically, we use many- 
body theory [5-7,11] (MBT) to calculate the annihilation 
7-ray spectra and annihilation rates for positron annihi¬ 
lation on individual subshells of the noble gases Ar, Kr 
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and Xe. We show that the theory is the first to accurately 
reproduce the the measured spectra [1] (see Ref. [10] for 
a brief account of this part of the work). Furthermore, 
we establish firmly the fraction annihilation on individual 
subshells, and uncover a simple scaling of the correlation 
enhancement factors, which describe the increase in the 
annihilation probability beyond the independent-particle 
approximation (IPA), with the electron ionization energy 
[10]. This work, taken together with our work in Refs. 
[5] and [10], gives a (near) complete description of the 
positron noble-gas system. 

The ability to discern the core signal in the 7-spectrum 
was first realised experimentally in [12]. The elemental 
specificity arising from the chemically distinct signature 
in the large Doppler-shift regions of the spectra corre¬ 
sponding to core annihilation has since been exploited 
[13-15], and measurement of the signal of individual se¬ 
lected core levels has been demonstrated [16]. Positron 
annihilation on core electrons is now a key process in sev¬ 
eral established experimental techniques, e.g., positron- 
induced Auger-electron spectroscopy (PAES) [17-20], 
used for diagnostics of industrially important materials, 
and the recently emerged time-resolved PAES [21], to 
study the dynamics of catalysis, corrosion, and surface 
alloying [22]. Coincident measurements of the 7-rays and 
Auger electrons yields 7-ray spectra for individual core 
orbitals [16, 23]. 

In all cases, interpretation of the measured annihila¬ 
tion 7-spectra relies heavily on theoretical input. For 
example, to interpret the results of PAES one needs to 
know the relative probabilities of positron annihilation 
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with inner electrons of various atoms [20]. These are of¬ 
ten computed using the independent-particle approxima¬ 
tion (IPA) [12, 24], or in a modified IPA framework that 
includes phenomenological enhancement factors that at¬ 
tempt to account for the important effects of electron- 
positron correlations [15, 16, 20, 23, 25-29]. In most 
cases, the materials of interest are condensed matter sys¬ 
tems, and the enhancement factors are usually calculated 
using a variety of density functional theory methods, e.g., 
the local density approximation (LDA) [2, 30], the gen¬ 
eralized gradient approximation (GGA) [31-34], or the 
weighted density approximation (WDA) [35-37], all of 
which rely heavily on theoretical considerations of the 
electron gas. However, due to the strong variations in 
the density, the LDA is not expected to work well for the 
core electrons, and has been found to overestimate the 
annihilation rates [37] The enhancement factor approach 
within the one-component and Boronski-Nieminen two- 
component LDA [30], as well as within the GGA, have 
been tested and compared with bound positron-atom 
stochastic variational calculations in Ref. [38], although 
no conclusive statements were made regarding the core 
enhancements. 

One particular weakness of the DFT enhancement fac¬ 
tor approach is that it obscures the true non-local many- 
body nature of the problem. The early theoretical works 
of Kahana, Carbotte and others (see, e.g., Refs. [39-49]) 
on positron annihilation were based on many-body the¬ 
ory, which inherently incorporates the non-local effects. 
Naturally however, these methods contained numerous 
approximations, and again were mainly annihilation rate 
calculations for the electron gas. Some effects of the non¬ 
local components of the amplitude have been included in 
phenomenological enhancement factor calculations using 
the WDA [37], and to a lesser extent in the GGA [15, 26]. 
However, these methods again rely on considerations of 
the electron gas. 

Regarding the core enhancement factors in particular, 
there are a number of calculations that differ in their 
treatment and overall, there is no definite agreement in 
their values. At one extreme it is assumed that the core 
enhancement is negligible, while at the other it is as¬ 
sumed to be equivalent to the valence enhancement (see, 
e.g., [37] and references therein). In reality, the core elec¬ 
trons are less susceptible to a perturbation from the inter¬ 
action with the incident positron than the valence ones. 
One should therefore expect the core spectra to be en¬ 
hanced by a factor 1 < 7 core <7vaiencej be., one should ex¬ 
pect the correlational enhancement to be weaker for the 
tightly-bound core than for the valence electrons. Indeed, 
a recent study by the authors that investigated the cor¬ 
relational enhancement for annihilation of positrons on 
H-like ions, for which the electron binding energies span 
a range typical of those of core electrons in many-electron 
atoms, suggested that the core enhancement could be sig¬ 
nificantly greater than unity [7]. 

For positron annihilation in condensed matter systems, 
an alternative approach to the calculation of core electron 


enhancement factors is suggested by the fact that the core 
electrons of these systems are localized and atomic in na¬ 
ture. Their enhancement factors are therefore open to ab 
initio atomic-based calculations. This was realised by 
Bonderup et al. [45], who, in order to elucidate the role 
of electron-positron correlations in condensed matter sys¬ 
tems, made a first-order perturbation theory estimate of 
the enhancement for positron annihilation on hydrogen¬ 
like ions. Compared to condensed matter systems, the 
problem of positrons annihilating in atomic systems is 
in principle, much simpler. For one, simplifications re¬ 
sult from the spherical symmetry of the atom that is not 
present in the condensed matter systems. For the latter, 
increased complexity comes from the periodic crystal po¬ 
tential: calculating the electronic and band structure is 
much more involved, and from the Bloch-wave nature 
of the incident positron wave function, which must be 
properly accounted for. Despite this, for annihilation 
on core electrons, positron-condensed-matter studies are 
much more numerous than corresponding positron-atonr 
studies. 

Evidence of the core contribution to the annihilation 
7 -spectra of the noble gases was first reported by Iwata 
et al.[ 1], in which upper bounds (2% and 3% for Ar and 
Kr respectively) were set on the core contribution to the 
total spectra. This was followed by the first experimental 
study of positrons annihilating with inner-shell electrons 
in isolated two-body positron-atom interactions [50], in 
which more precise measurements of the 7 -spectra of 
the noble gases Ar, Kr and Xe were reported. In that 
work, a static Hartree-Fock calculation was successfully 
used to identify the inner-shell electron contribution to 
the measured spectra. However, owing to the fact that 
it ignored the state dependent enhancement, the rela¬ 
tive fraction of annihilations with inner-shell electrons 
was not predicted accurately by this theory. Almost a 
decade later, Dunlop and Gribakin [ 6 ] used many-body 
theory to explain the main features of the 7 -spectra of 
the noble gases. In that work, however, only annihila¬ 
tion on the valence shells was considered, the short-range 
electron-positron interaction was treated only in first- 
order, and positron-atonr correlations were ignored; sig¬ 
nificant discrepancies between the theoretical results and 
experiment were evident. Dunlop [51] later attempted 
to calculate the core contribution to the spectra via a 
semi-empirical fitting procedure. For this, the annihila¬ 
tion on the valence shells was calculated with the full in¬ 
corporation of short-range electron-positron correlations, 
whereas the core was treated in the simple independent- 
particle-model (IPA) approximation. It will be shown 
however, that neglecting the core enhancement in this 
way leads to an overestimate of the core contribution to 
the 7 -spectrum. 

The purpose of this work is three-fold: (1) to examine 
in detail the effect of electron-positron correlations on the 
process of annihilation on tightly-bound core electrons of 
many-electron atoms (using the noble gases as test cases) 
and in doing so; ( 2 ) fully explain the measured spectra of 
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the noble gases and establish the fraction of annihilation 
on individual core and valence subshells, and finally; (3) 
to calculate the ‘ true ’ core-electron enhancement factors 
(and demonstrate that they can be used to correct rela¬ 
tively simple independent-particle-approximation calcu¬ 
lations to produce accurate spectra). 

The structure of the paper is as follows. In Sec. II we 
outline the many-body theory of annihilation 7 -spectra, 
detailing the form of the fully correlated incident positron 
wave function and annihilation vertex. In Sec. Ill we 
briefly discuss the numerical implementation of the the¬ 
ory. In Sec. IV we present our calculated annihilation 
spectra for Ar, Kr and Xe, comparing them with the 
experimental results of Iwata ef al. [50]. In addition, 
the partial Z c g are calculated for each subshell. In 
the final part, we calculate vertex enhancement factors 
and demonstrate how they can be used in independent- 
particle-model calculations to obtain accurate annihila¬ 
tion spectra. A summary is given in Sec. V. 


II. MANY-BODY THEORY OF ANNIHILATION 
7-SPECTRA 

A. Basics 

The process of electron-positron annihilation is de¬ 
scribed fundamentally by quantum electrodynamics [52, 
53]. It dictates that low-energy positron annihilation 
proceeds predominately via two photon production, a 
process in which the total spin of the electron-positron 
system must be zero [52, 53]. Consider then a positron 
with momentum k annihilating in a many-electron sys¬ 
tem with an electron in state n [54] . In the annihila¬ 
tion event, two gamma-photons with total momentum 
p 7l + p 72 = P are produced. Let us denote the ampli¬ 
tude for this process by 

A„ k (P) = (T^- 1 ; 27 P|^ +1 ), (1) 

where 4/^ r+ 1 (ri,... rjv; r) is the true, fully-correlated 
many-particle wave function of the initial state of N elec¬ 
trons (coordinates 17 ) and the positron (coordinate r), 
and T() r_ 1 (r 1 ,... rjv-i) is that of the final state of sys¬ 
tem with N — 1 electrons and a hole in state n. At large 
positron-atom separations the initial N + 1 particle wave 
function has the form 

< + 1 (r 1 ,...r JV ;r)^<(r 1 ,...r JV )e ikr ( 2 ) 


is Doppler shifted. For example, for the first photon 
I? 7l = E 1 + mcV cos 9, where V = P/2 m is the centre- 
of-mass velocity of the electron-positron pair and 9 is the 
angle between the direction of the photon velocity c and 
V. Assuming that V -C c, and p 7l = E 7l /c « me, the 
shift of the photon energy from the centre of the line, 
e = £? 7l — E^, is 

Pc 1 

e = me V cos 9 = — cos 9 = P • c. (3) 

2 2 

The typical momenta of electrons bound with energy e n 
determine the Doppler width of the annihilation spec¬ 
trum, e ~ Pc ~ y/\e n \mc 2 |e n |. Hence the shift of the 
line centre e n /2 from E~ = me 2 = 511keV can usually 
be neglected, even for the core electrons. The annihila¬ 
tion 7 -ray (Doppler) spectrum can then be written in a 
form similar to a Compton profile, as the projection of 
the annihilation probability: 

m nk (e,c)= J K k (P)| 2 «5^-ip-c) (4) 

For gaseous systems, however, the experimentally acces¬ 
sible quantity is the average of this over the direction of 
emission of the annihilation photons c = c|c|. Also, to 
compare with experiments that use positrons confined in 
a trap, one must also average over the direction of the 
incident positron momentum k giving 


w ne (e) 



c ) 


dfl c dff k 
47T 47T 


1 r°° r 

C J 2|e|/c JS 2 P 


K e (p )| 2 


dflp 

( 2 ^ 


PdP 


( 5 ) 
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where |A r!£ (P )| 2 = J |A„ k (P)| 2 ciH k / 47 r. The annihila¬ 
tion amplitude A(P) is therefore central to the evalua¬ 
tion of the annihilation 7 -spectra. Its calculation from 
many-body theory is outlined below. In most experi¬ 
ments, including those to which we compare our results, 
the annihilation photons are not detected in coincidence 
with the final state of the atom. The observed spectrum 
is then given by a sum over the electronic states of the 
system: w e (e) = Y n ™ ne ( e ) m As we shall see, the total 
spectrum still retains distinct characteristics of the indi¬ 
vidual final states. 

A useful quantity that is related to the 7 -spectra is the 
normalized annihilation rate parameter Z e g. It is defined 
as the ratio of the true annihilation rate A to the spin- 
averaged Dirac annihilation rate on a free electron gas of 
number density n e , Xd = ^r^cne [55, 56] 


which is normalized to the positron plane wave, where 
'S’q is the atomic ground state wave function. In the 
centre-of-mass frame, where the total momentum P is 
zero, the two photons propagate in opposite directions 
and have equal energies E y = p-yC = me 2 + \ {E\ ~ Ef) ~ 
me 2 ~ 511 keV, where E-, and Ef denote the energy of 
the initial and final states (excluding rest mass). When 
P is non-zero however, the two photons no longer prop¬ 
agate in exactly opposite directions and their energy 


A = X D Z eS , (7) 

where ro = e 2 /mc 2 is the classical radius of the electron 
(in CGS units). If the electrons belong to the atoms 
in a gas of number density n, then n e = Zn, where Z 
is the number of electrons per ion. The corresponding 
spin-averaged annihilation cross section is 

- _ A _ 2 c 7 

7 — — TT^O ^eff • 

n e v v 


( 8 ) 
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where v is the positron velocity. By definition, Z e g gives 
the effective number of electrons that contribute to the 
annihilation. In this way it provides a measure of the 
strength of the correlations. It can be determined from 
the annihilation 7 -ray Doppler spectrum through the re¬ 
lation 

/ OO 

w(e)de. (9) 

-OO 

The measurements to which we will compare our calcula¬ 
tions are those of Surko and co-workers of the San Diego 
group [50]. These experiments used a positron trap and 
low gas densities. In this regime, the parameter Z e g is 
independent of n e . 

B. Many-body expansion of the annihilation 
amplitude 

The many-body theory of positron annihilation 7 - 
spectra was introduced in [ 6 ]. In that work however, 
only the zeroth- and first-order correction to the anni¬ 
hilation vertex were discussed and the important effect 
of positron-atom correlations on the incoming positron 
wave function were neglected. In this section we give a 
concise overview of the theory. We present formulae for 
the full annihilation vertex and briefly discuss the treat¬ 
ment of the modified positron wave function. 

Taking the non-relativistic limit of the two-photon 
QED annihilation amplitude leads to the second quan¬ 
tized operator of two-photon annihilation O a ( P) [ 6 , 57- 
59] 

Oa(P) = J d 3 r e-* p ->(r)^(r), (10) 

where ^(r) and f>(r) are the electron and positron field 
annihilation operators respectively [60] . In this form, it 
is clear that non-relativistic annihilation occurs when the 
positions of the particles coincide. To apply the many- 
body theory, a complete orthonormal set of single particle 
wave functions must be introduced. The expansion takes 
its simplest form in Hartree-Fock (HF) single-particle ba¬ 
sis | fi), in which the held operators of the electron and 
positron become 

= ; £( r ) = X^( r ) ( n ) 

n n 

where ip^fr) and <^(r) are the HF wave functions (of 
central-held type), and a^ and b^ are the HF held anni¬ 
hilation operators, of the electron and positron. In this 
basis the operator of annihilation becomes: 

O a (P) =^2(P\5\^u) a^b v , (12) 

where we use the notation 

(P| S\nv) = J e~ lP - r iJjffr)ip u (r)d 3 r (13) 


The Hamiltonian of the system is 

H = H HF + H intl (14) 

where Huf is the Hartree-Fock Hamiltonian that satishes 
Hhf\^o) = E o F \®o), where \®o) = a=i«l| 0 ) is the 
IV-electron ground-state HF determinant wave function. 
The many-body expansion of A ne ( P) proceeds through 
the interaction Hamiltonian H mt = I4es + O a ( P), where 
the residual interaction V les = V e - e + + V e - e - — U^ F is 
the sum of the standard two-body electron-positron and 
electron-electron Coulomb operators minus the one-body 
HF operator (for the frozen core of the iV-electron sys¬ 
tem). 

The resultant diagrammatic expansion is shown in 
Fig. 1 (see the figure caption for an explanation of the di¬ 
agrams). The construction of the diagrams follows stan¬ 
dard rules. First, note that the use of HF single-particle 
basis ensures that only pure Coulomb operators appear 
in the final diagrams, i.e., there are no explicit contri¬ 
butions from Chf- Furthermore, since there is a single 
incident positron, only terms in the expansion involving 
a single O a (P) are non-zero [61]. Single-hole final-state 
diagrams at order m = nil + m 2 in the Coulomb inter¬ 
action (wavy line), where mi is the number of electron- 
positron Coulomb interactions and m 2 is the number of 
electron-electron Coulomb interactions, have mi positron 
propagators (single directed lines) and mi + 2 m 2 elec¬ 
tron propagators, in addition to a single positron and 
hole wave function (labelled e and n respectively), and 
the two 7 -rays (double-dashed lines). 


1. The annihilation vertex 

Figure 2 shows the main contributions to the annihi¬ 
lation vertex: the zeroth-order (IPA) vertex, as well as 
the non-local first-order and the positron-electron ladder 
series, or T-block’, corrections. From inspection of the 
diagrammatic expansion, one can write the exact annihi¬ 
lation amplitude in the abstract form 

A„ e (P) = [ e~ lPr {if e (r)ip n {r) 

{ ( 15 ) 
+ Ap(r;ri,r 2 )^ £ (ri)</J„(r 2 )dridr2 j dr, 

where the first term is the zeroth-order (IPA) vertex 
(Fig. 2a), and Ap is the non-local annihilation kernel. 
In this form, the amplitude A(P) is clearly seen to be 
a Fourier transform of the correlated pair wave function 
(the term in the brace), whose modulus squared can be 
referred to as the electron-positron annihilation momen¬ 
tum density. In the HF basis it takes the explicit form 
(using the notation of Eq. 13) 

A„ E (P) = (P|(g)|n, £ ), 


(16) 
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FIG. 1. Diagrammatic expansion of the annihilation amplitude A ne (P). Diagram (a) represents the zeroth-order (independent- 
particle-model) approximation. All other diagrams are corrections to this that can be grouped into two independent types. The 
first type, which includes diagrams (b), (c) and (f) are corrections to the annihilation vertex. The second type, which includes 
diagrams (d) and (e) are corrections to the zeroth-order Hartree-Foc.k wave functions. The general form of the amplitude, shown 
on the LHS of the equality, consists of four elements: (1) the incident positron wave function labelled by the quantum numbers 
e (the double line indicates the ‘dressed’, i.e., fully correlated wave function); ( 2 ) the general annihilation vertex (crossed circle); 
(3) the wave function of the annihilated electron (hole) of state n, and finally; (4) the two 7 -rays of total momentum P. 




FIG. 2. Main vertex contributions to the annihilation amplitude A ne ( P). Shown on the right hand side are the different orders 
of the annihilation vertex: (a) zeroth-order (IPA); (b) first-order correction, and; (c) T-block’ correction that contains the 
electron-positron ladder series. 


where the vertex operator 


in Fig. 2) 


= *-£ 


fi,iy 


S\n,v){v,lAV 

£ + £n — £u + £ 1 / 


£ 


%2,^2)(^2,M2|r £+£ jMl,ZA)(t / l 1 Ml|T / 

(e + e n — £fj , 2 — £i/ 2 )(e + £ n ~ £/ii — £1/1) ’ 


(17) 


where the sum over /i and v is over the excited electron 
and positron states. Equation 16 is written to emphasize 
the components common to all diagrams and is in accor¬ 
dance with the general diagrammatic form shown in Figs. 
1 and 2 : the external positron wave function (labelled e), 
linked through the general vertex 0 to the annihilated 
electron (hole) wave function (labelled n) and the two- 
gamma rays P. Note that the zeroth-order (IPA) result is 

just that of Eq. 13, i.e., An] (P) = (P|<j|ne). The T-block 
is determined from the linear integral equation, arrived at 
trivially by inspection of its diagrammatic form (shown 


(a/3|r B |7<5) 


(aj3\V\^5) + 


((3\V\nv)(vti\T E \'y8} 

E - - e v + ir] 


(18) 


where r) is a number close to zero. The successful com¬ 
putation of the T-block correction to the vertex was first 
performed in this way in [11], in which calculations of Z e g 
for positron annihilation on hydrogen were reported. It 
was later incorporated in the calculations of the 7 -spectra 
for annihilation on valence shells of the noble gases [51], 
and recently to annihilation on H-like ions [7]. It can 
give rise to a significant enhancement above that of the 
first-order correction. 

The Hartree-Fock wave functions involved in the 
matrix elements are of central field type (fintmi r ) = 
}.P n t:(r)Y( rn (r). The incident positron wave function is 
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a continuum state of the form 

V’k(r) = ^^J^J 2 i eelSe Y; m (k)Y em (r)P ee (r), 

(19) 

where P(r ) is normalized to a delta function of energy 
(in Rydbergs). The positron radial wave function P(r) 
is modified from the HF one by the positron-atom corre¬ 
lations, as explained below. 


2. Positron ‘Dyson’ quasiparticle wave function 

The correlated incident positron quasiparticle wave 
function ijj e can be calculated from the integral form 
of the Dyson equation for the positron Green’s function 
(see, e.g., [62, 63]): 

(h 0 + E £ ) ip e (r ) = eip £ (r) , ( 20 ) 

where Hq is the unperturbed Hamiltonian of the N + 1 
particle system, which we take in this work to be the 
Hartree-Fock Hamiltonian, and E £ is the non-local cor¬ 
relation potential. It is an integral operator that acts as 
S £ ^ e (r) = / S £ (r, r')ip e (r')dr' , and is equivalent to the 
energy-dependent irreducible self-energy of the positron 
in the presence of the atom. The main contributions 


e v e' £ e' £ £' e 



FIG. 3. The main contributions to the irreducible positron 
self-energy E, or more specifically, the matrix elements 
(e|E £ | e '>. 


to E are shown in Fig. 3. At large positron-atom dis¬ 
tances E reduces to the local dipole polarization potential 
E « —T^S(r — r'), where ad is the dipole polarizability 
of the atom in the HF approximation 



n,/x 


M I »I 2 


( 21 ) 


Formally, the quasiparticle wave function (or ‘Dyson 
orbital ’) if> e is given by 


V’e(r) = f % N (ri, • ■ ■, r A r)«'^ + 1 (ri,..., r N ; r)d N r 1 

( 22 ) 

i.e., it is the projection of the fully-correlated N + 1 
particle wave function 'F JV+1 (of energy E) onto the N- 
electron atomic ground state of N (of energy Eq), 


and e = E — Eq. It can be calculated as outlined in Ref. 

[ 11 ], as 

V>e(r) = ¥ , i° ) ( r ) + P j <M r ) de '■ ( 23 ) 

where E is the reducible self-energy and where ‘P’ means 
the principal value of the integral. To normalise this 
wave function to a delta-function in energy (Rydbergs), 
it must be multiplied by the normalisation factor 

N(e)= (l + ( £ |E £ | £ ) 2 )" 1/2 . (24) 

Overall, this method of determining the positron quasi¬ 
particle wave function is equivalent to solving the 
Schrodinger equation with an optical potential [64]. 


III. NUMERICAL IMPLEMENTATION 

The numerical procedure for evaluating the dia¬ 
grams has been successfully demonstrated previously in 
Refs. [11] and [ 6 ], where it was used to calculate Z e g for 
annihilation on atomic hydrogen, and the 7 -spectra for 
annihilation on the valence shells of the noble gases. The 
electrons involved in these studies have relatively small 
ionization energies, compared to the core-electrons that 
are of interest here. More recently, the authors performed 
a study of positron annihilation on hydrogen-like ions [7] . 
The ions considered there had electron ionization energies 
similar to those of typical core electrons of many-electron 
atoms. Using other, variational results as a benchmark, 
this study provided verification of the applicability and 
accuracy of the numerical method to annihilation with 
electrons of relatively large binding energies. Full details 
of the numerical procedure can be found in the affore- 
mentioned works. Here, we aim to reduce repetition to 
the minimum compatible with the purpose of this section. 

The main quantities to be calculated are the self¬ 
energy, the annihilation amplitude, and the 7 -spectra 
w(e), whose angular reduced formulae are given in Ap¬ 
pendix A and Ref. [11]. To evaluate the various many- 
body diagrams, one must first generate a set of single¬ 
particle Hartree-Fock basis states. For the electron, this 
basis consists of the atomic ground-state wave functions 
as well as excited electrons in the field of the A^-electron 
frozen core of the atom. For the positron, the static- 
field positron-atom potential is repulsive, and all states 
are continuum states in the field of the A-electron frozen 
core of the atom. To evaluate these, we use a mesh in 
momentum space of 201 points starting from threshold 
with a step size of k = 0.02 a.u. 

To facilitate the calculation of the diagrams, partic¬ 
ularly those involving the ladder block, we introduce a 
B-spline basis [ 6 , 7, 11, 65-67]. B-splines of order k are 
piecewise polynomials of degree k — 1 defined over a re¬ 
stricted domain (‘box’) that is divided into n — k + 1 
segments by a knot sequence of n points: r,; £ [0, R], 
where n is the number of splines in the basis [65] . In this 
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work, we use an exponential knot sequence defined as 
ri = po{e al — 1), where po = 10 -3 a.u. and a is determined 
by the condition r n -k+i = R, the box size. It allows for 
the accurate description of both the ground-state atomic 
wave functions, which can have many oscillations over 
the typical size of the atom and rapidly vanish outside 
the atom, as well as the continuum states that extend 
to larger distances up to the confinement (box) radius. 
The choice of the exponential knot sequence also ensures 
rapid convergence in the summations over intermediate 
states in the diagrams. The choice of spline basis param¬ 
eters and box size is discussed below. 

By expanding the single-particle HF orbitals P n t(r) = 
JT c"^-Bj(r), where are the expansion coefficients, 
and Bi is the i- th spline of the basis, the Schrodinger 
equation can be reduced to the discrete eigenvalue prob¬ 
lem HWc e = eBc e , where ‘ B i {r)H^' 1 Bj{r)dr 

and Bij = Bi(r)Bj(r)dr. Note that to implement 
the boundary conditions P( 0) = P(R) = 0 the first and 
last splines are discarded in the expansions. The solu¬ 
tions of the equation for a given angular momentum i 
are a set of n — 2 eigenfunctions. For the electron, the 
lowest energy states from this set correspond directly to 
the ground-state wave functions. The remaining electron 
states of the set describe excited electrons propagating in 
the frozen-core Hartree-Fock field of the N ground-state 
electrons. The solutions for the positron are continuum 
states that describe the positron propagating in the field 
of the HF frozen core of the iV-electrons. 

The choice of optimal spline basis parameters and box 
size depends on a number of competing considerations. 
These have been discussed in detail in Refs. [7] and [11]. 
For example, the vertex corrections to the annihilation 
amplitude involve small electron-positron distances and 
their accurate evaluation therefore requires an adequate 
minimum spatial resolution. Increased resolution can be 
obtained by either using a basis with a larger number of 
splines, or by reducing the box size, or both. For a given 
box size, however, increasing the number of splines in the 
basis reduces the density of states: a greater number of 
states must be included in the numerical calculations to 
achieve the same energy coverage and reach convergence. 
On the other hand, the positron-atom correlations in¬ 
volve are of longer-range and their accurate evaluation 
requires a large box size. 

The calculations were performed using a basis of n = 
40 splines of order k = 6, and repeated using a basis of 
n = 60 splines of order k = 9. For both spline bases, 
two different box sizes, R = 20 a.u. and R = 30 a.u., were 
used. We found that the n = 60, k = 9 core-subshell 
results suffered from poor convergence, and that the best 
overall results were obtained using the basis of n = 40 
splines of order k = 6 with a box size of R = 20 a.u. 
The lowest 20 B-splines of the basis were included in the 
summations over intermediate states, corresponding to 
an energy coverage from zero to > 10 2 a.u., significantly 
greater than the energy of the valence and subvalence 
(core) ground-state wave functions that we are interested 



FIG. 4. Absolute values of the energy of B-spline basis states 
(n = 40 splines of order k = 6) with angular momentum 
1 = 0 for electron (empty circles) and positron (solid circles) 
for Ar, Kr and Xe. The l = 0 electron ground states are 
labelled. For Ar, the energies of the s-states are 118.61, 12.32 
and 1.28 a.u. For Kr, the energies of the s-states are 520.17, 
69.90, 10.85 and 1.15 a.u. For Xe, the energies of the s-states 
are 1224.40, 189.34 and 40.18 a.u. The ground p- and d-states 
are also shown as labelled squares and diamonds, respectively. 
Spline states with l > 0 have larger energies than those with 
l = 0 . 





in (see Fig. 4). All results quoted henceforth have been 
calculated with these parameters. 

When evaluating matrix elements involving the inci¬ 
dent positron wave function e, e.g., the self-energy ma¬ 
trix (e|E e |e), or the annihilation amplitude A„ e (P) = 
(P|0|n,e), we use the true HF continuum states |e). 
This can be done by making further use of the spline 
completeness, introducing a resolution of the identity to 
re-write the matrix elements as, e.g., 

(e'lE^le) = J2{e'\f\v){v\f- 1 T, E f- 1 \i'')(i''\f\£),(25) 

is,is' 

(P|%,e) = ^(P|<5/- 1 |z/)(z/|/|n,e). (26) 

V 

The insertions of / _1 /, where / = (R — r), are made to 
nullify any numerical error arising due to the fact that 
the B-splines are zero at the boundary R, whereas the 
Hartree-Fock states are not. In this way, the matrix ele¬ 
ments of all quantities involving the true positron states 
can be evaluated as matrix elements involving individual 
spline basis states. 

In the calculation of the positron self-energy (Fig. 3) we 
include only those hole states n in the diagrams that cor¬ 
respond to the two valence electrons, owing to the negligi¬ 
ble contribution to the dipole polarizability of the atom 
ad from the strongly-bound core electrons. Any small 
effect of omitting the core states from the self-energy di¬ 
agrams will manifest only in the positron wave function, 
and should have a negligible effect on the vertex enhance¬ 
ment. In the diagrams of the annihilation amplitude the 
final hole state n is, of course, taken to be any of the 
ground states, including the core subshells. 

The finite box size leads to a cut-off in the self-energy 
matrix element: it is evaluated only over the domain of 
the box. To compensate for this, we make use of the 
long-range asymptotic form of the correlation potential 
and add to the spline-based matrix element (calculated 
over the domain of the box) the contribution at r > R: 

J r P sA r ) (-^f) P e'iA)dr, (27) 

where the asymptotic radial wave functions P e are calcu¬ 
lated with the HF phase shift as 

p e i e {r) ^ [ji e (kr ) cos^ F - n h (kr ) sin^ F ] ,(28) 

where ji and ni are the spherical Bessel and spherical 
Neuman functions, respectively. 

For the calculation of the annihilation momentum dis¬ 
tribution we used a mesh in the y-ray momentum of 101 
points from 0 to 9a.u. We also repeated the calcula¬ 
tions with 101 points from 0 to 15a.u., although no dis¬ 
cernible difference was found in the resulting y-spectra 
in the range of Doppler-shifts of interest (0-9 keV). 

Finally, convergence with respect to the highest an¬ 
gular momentum of the intermediate excited states used 


in the calculation, f max —> oo, was achieved using the 
asymptotic formula [6, 7, 68] 

(e) = (e) + ^max + l/2)-\ (29) 

where, for a given k, A is a constant that is determined 
from the results of the calculations performed using dif¬ 
ferent f max . In this work all diagrams were calculated for 
4nax = 7, 8, 9 and 10. 


IV. RESULTS AND DISCUSSION 

In this section we present the results of the MBT cal¬ 
culations. In Sec. IV A we discuss the effect of the corre¬ 
lations on the spectra of individual core and valence sub¬ 
shells of the noble gases Ar, Kr and Xe, and compare our 
calculated annihilation y-spectra with the existing exper¬ 
imental spectra. In addition, we calculate the total, and 
fractional contributions to, the normalized annihilation 
rate parameter Z b r for each subshell. In Sec. IV B, we 
use the results of the ab initio many-body theory calcu¬ 
lations to provide state-dependent enhancement factors 
for the core and valence electrons that quantify the vertex 
enhancement to the annihilation vertex. All the results 
presented have been calculated using an s-wave incident 
positron of thermal momentum (k = 0.04a.u.). 


A. Annihilation y-spectra of noble gases 

1. Many-body calculated annihilation 'y-spectra of valence 
and subvalence shells of Ar, Kr and Xe 

The effects of the correlations on the y-spectra are 
shown in Fig. 5, for the specific examples of the valence 
and subvalence p-subshells of Ar, Kr and Xe, respectively. 
For each subshell, the spectra calculated using both the 
HF (static) and Dyson incident positron in the three 
different approximations of the annihilation vertex (see 
Fig. 2) are shown. As the spectra are symmetric about 
zero energy shift, only positive Doppler energy shifts are 
shown. Interpretation of the spectra is aided by Table I, 
which gives the values of Z e g calculated as an integral of 
the spectra through Eq. 9. 

For each atom, for a given annihilation vertex order, 
the spectra calculated using the Dyson positron wave 
function are significantly greater in magnitude than those 
calculated using the HF one, e.g., the spectrum of the 
valence np subshell of Ar is increased above its Hartree- 
Fock result by ss 8 times; for Kr by ss 18 times, and for 
Xe by as much as ss 70 times. The reason for this wave 
function enhancement is the introduction of an attrac¬ 
tive potential through the strong positron-atom correla¬ 
tions. This potential supports low-lying positron virtual 
s levels, which give rise to the characteristic resonant be¬ 
haviour: Z e s(k) oc (k 2 + A; 2 ) -1 , where k = 1/a is the 
reciprocal of the scattering length a, assumed here to 


y-spectra (keV" ) y-spectra (keV" ) y-spectra (keV" ) 
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FIG. 5. y-spectra for valence and subvalence p-subshells of Ar, 
Kr and Xe, calculated using HF (static) and Dyson positron 
wave function in different orders of the annihilation vertex. 


be much larger than the radius at which annihilation oc¬ 
curs (typically the size of the atom) [69-72]. The relative 
strength of the wave function enhancement of Ar, Kr and 
Xe corresponds to the increase in their respective scatter¬ 
ing lengths: for Ar a ~ — 4.4a.u., for Kr a ~ —lO.la.u., 
and for Xe a ~ —81 a.u. [73]. Incidentally, this resonant 
behaviour means that Z e g drops off dramatically as k 
increases, and thus a proper comparison with the exper¬ 
imental Z e s measured using thermalized positrons can 
only be made after performing a Maxwellian average of 
the calculated Z c g(k) over k. Since it is the calculation 
of the y-spectra that is the main interest here, this fact 
is not of utmost concern. Nevertheless, for comparison’s 
sake the Maxwellian averaged Z e g of the valence shells, 
calculated previously using MBT [73], are quoted in the 
table. 

For a given atom, deeper lying subshells experience a 
greater wave function enhancement. For example, for a 
given vertex order the Dyson spectrum for the Ar valence 
3 p subshell is « 8 times larger than the HF spectrum, 
compared to the subvalence (n — l)s, for which it is ss 
13 times larger. For Xe, the valence 4 p subshell is ss 
70 times larger than the HF spectrum, compared to the 
subvalence (n — l)s, for which it is ~ 160 times larger. 
This is a result of the fact that the positron must tunnel 
through a repulsive potential to reach the core electrons. 
The additional attraction from the positron-atom corre¬ 
lations effectively lowers the HF repulsive potential: the 
wave function dies away more slowly in the tunnelling re¬ 
gion, and consequently, the overlap of the positron wave 
function and the core electrons is increased. Compara¬ 
tively, the overlap of the positron and valence electron 
wave functions is relatively unchanged. 

Concentrating now on the effect of the vertex correc¬ 
tions, it is clear that the addition of the first and F- 
block vertex corrections enhances the spectra above the 
zeroth-order IPA result. Fundamentally, this is due to 
the increase of the annihilation probability owing to the 
stronger attraction between the annihilating pair. This 
attraction is of particular importance, given the fact that 
non-relativistic annihilation occurs when the positions of 
the positron and electron coincide. For both the core 
and valence shells, the addition of the first-order cor¬ 
rection to the IPA vertex increases the magnitude of 
the spectra considerably: e.g., for the 4 p valence and 
3 p subvalence subshells of Kr the spectra is increased 
at zero energy shift respectively by ~ 2.5 times from 
4.36 to 11.44 keV -1 , and by ~ 1.3 times from 0.025 to 
0.034 keV -1 . For the valence shells, the addition of the 
higher order F-block correction causes a further signif¬ 
icant enhancement: e.g., for the 4 p subshell of Kr, it 
produces a further 2.3 times (from 11.44 to 27.6 keV -1 ) 
enhancement to the first-order result. For the core how¬ 
ever, the result of adding the T-block correction is much 
less dramatic, only raising the spectra slightly above that 
obtained with the first-order correction: e.g., for the 3 p 
subshell of Kr, the T-block correction only increases the 
magnitude of the first-order result by a further 1.05 times 
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TABLE I. Calculated Z e g (using Eq. 9) and 7 -spectra FWHM for annihilation on the valence n and subvalence (n— 1) subshells 
of the noble gases. Values are given for the different approxmations to the incident positron: HF orbital and Dyson quasiparticle 
wave function; and to the annihilation vertex: 0th, 0+lst and 0+lst+T (Fig. 2). The values are for an s-wave incident positron 
of momenta k = 0.04 a.u. 


Subshell 



Z e ft 


7 -spectra FWHM (keV) 



HF 


Dyson 

HF 

Dyson 


0 

0+1 

0+1+F 0 

0+1 0+1+r 

0 0+1 0+1+r 

0 0+1 0+1+r 


Argon 


np 

0.606 

1.530 

3.144 

5.651 

13.43 

24.97 

2.88 

2.71 

2.40 

3.12 

2.95 

2.63 

ns 

0.136 

0.260 

0.343 

1.387 

2.576 

3.320 

1.87 

1.82 

1.77 

1.98 

1.93 

1.88 

(n — 1 )d 













(n - 1 )p 

0.622 [- 2 ] 

0.849[-2] 

0.891 [-2] 

0.812[-1] 

0.110 

0.115 

9.42 

9.28 

9.24 

9.57 

9.43 

9.39 

(n — l)s 

0.235[-2] 

0.307[-2] 

0.319[-2] 

0.308[-l] 

0.400[-l] 

0.414[-1] 

5.19 

5.20 

5.20 

5.27 

5.27 

5.27 

Total 

0.75 

1.80 

3.50 

7.15 

16.16 

28.44 






2.52 

(n tot.)fe a 






26.5 







Expt. b,c 






26.77+0.09 

Krypton 






2.30 

np 

0.565 

1.543 

3.742 

13.12 

16.48 

68.67 

2.56 

2.41 

2.12 

2.81 

2.67 

2.37 

ns 

0.115 

0.231 

0.320 

3.018 

2.898 

7.744 

1.67 

1.63 

1.57 

1.78 

1.75 

1.70 

(n — 1 )d 

0.155[-1] 

0.233[-l] 

0.256[-l] 

0.510 

0.758 

0.827 

8.85 

8.70 

8.64 

9.09 

8.95 

8.88 

(n - l)p 

0.599 [-2] 

0.805 [-2] 

0.843[-2] 

0.205 

0.271 

0.283 

7.65 

7.58 

7.56 

7.77 

7.70 

7.68 

(n — l)s 

0.172 [- 2 ] 

0.221 [- 2 ] 

0.229[-2] 

0.591[-1] 

0.754[-l] 

0.779[-l] 

4.35 

4.35 

4.35 

4.42 

4.41 

4.41 

Total 

0.703 

1.81 

4.10 

16.91 

20.48 

77.6 






2.89 

(n tot.)*, a 






66.4 







Expt. b,c 






64.6+0.08 

Xenon 






2.10 

np 

0.510 

1.535 

4.725 

50.806 

136.67 

334.01 

2.23 

2.09 

1.81 

2.48 

2.37 

2.07 

ns 

0.104 

0.228 

0.351 

11.98 

24.64 

35.77 

1.41 

1.38 

1.33 

1.55 

1.51 

1.46 

(n — 1 )d 

0.229[-l] 

0.387[-l] 

0.450[-l] 

3.525 

5.671 

6.602 

6.73 

6.65 

6.60 

6.96 

6.89 

6.84 

(n - l)p 

0.628 [- 2 ] 

0.892[-2] 

0.952[-2] 

1.010 

1.403 

1.490 

5.79 

5.77 

5.76 

5.91 

5.88 

5.87 

(n — l)s 

0.156[-2] 

0 . 210 [- 2 ] 

0 . 220 [- 2 ] 

0.253 

0.335 

0.350 

3.35 

3.36 

3.36 

3.41 

3.41 

3.41 

Total 

0.645 

1.81 

5.13 

67.57 

168.0 

378.2 






2.01 

(n tot.}*, a 






402 







Expt. d,c 






401+20 






1.93 


a Thermally averaged valence total from many-body theory, including additional diagrams [73]. 
b Z e g total measured with thermalized positrons, from Ref. [74]. 

c Experimental 7 -spectra FWHM from single-gaussian fit of Ref. [1], Note that a more accurate description of the experimental spectra 
is given by the two-gaussian fit of Ref. [1], 
d Z e ff total measured with thermalized positrons, from Ref. [75]. 


(from 0.0340 to 0.0356 keV -1 ). A similar dominance of 
the first-order correction over the T-block correction was 
also seen for annihilation with the tightly-bound elec¬ 
trons of the H-like ions in [7]. Physically, the T-block de¬ 
scribes virtual positronium formation by the tunnelling of 
an electron from the atom to the positron. This process 
is more difficult for the tightly bound spatially compact 
core than for the valence electrons. 

From the figures it is also clear that for a given sub¬ 
shell, the enhancement from the vertex corrections is 
similar whether the HF or Dyson positron wave func¬ 


tion is used. Physically, this is to be expected, as the 
vertex corrections are to a large extent independent of 
the positron-atom correlation. We will discuss the ver¬ 
tex enhancement further in Sec. IV B, where we evaluate 
enhancement factors that quantify the important vertex 
enhancement and show how they can be used in simple 
IPA type calculations to calculate accurate annihilation 
7 -spectra. 

The correlations, both wave function and vertex type, 
evidently have a dramatic effect on the magnitude of 
the 7 -spectra and the annihilation rate. Their effect 
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FIG. 6 . The effect of electron-positron correlations on the 7 - 
spectrum FWHM for the 5 p subshell of Xe. Calculated using 
Dyson (heavy lines) and HF (light lines) positron in different 
orders of the annihilation vertex: zeroth (dashed); zeroth+ 
first (dash-dot line) and full (solid line). 


on the shape of the spectrum is less so, although it is 
still significant. Fig. 6 shows the spectrum of the Xe 
5 p subshell in the different approximations to the ver¬ 
tex and positron wave function. For ease of compari¬ 
son, all curves in the figure have been scaled to unity at 
zero energy shift. In addition, values of the full-width 
at half-maximum for all of the individual subshells cal¬ 
culated in the different approximations are given in Ta¬ 
ble I. For the valence subshells, it is clear that increas¬ 
ing the self-energy corrections in the wave function for 
a given vertex order affects a broadening in the spec¬ 
tra: the attractive positron-atom correlations accelerate 
the positron towards the electron resulting in a relatively 
large momentum distribution, and hence a broader 7 - 
spectrum. In contrast, increasing the vertex order affects 
a narrowing of the spectra: the first- and higher-order 
annihilation vertex corrections involve excited state (vir¬ 
tual) electrons that are relatively diffuse and therefore 
produce a relatively narrower momentum distribution, 
and hence a narrower 7 -spectrum. The degree of can¬ 
cellation that occurs between this competing narrowing 
and broadening leads to a relatively small change in the 
full-width at half-maximum (FWHM) between the sim¬ 
plest zeroth-order-vertex HF-wave function result and 
the full MBT (0-t-l+T)-vertex Dyson-wave function re¬ 
sult. This explains why one can obtain reasonable frac¬ 
tional core contributions through semi-empirical fitting 
of the zeroth-order HF-wave function calculated core and 
valence spectra to the total [50]. In contrast, the FWHM 
of the deeper lying subshells, e.g., the (n — l)s subshell, 
does not change noticeably. In summary, overall there is 
clearly some dependence of the FWHM on the approxi¬ 
mation used, and reliable core contributions to the total 


spectra can therefore only be calculated by taking full 
account of the correlations. 

Fig. 7 shows the 7 -spectra for the valence n, subva¬ 
lence ( n — 1) and (n — 2) shells for Ar, Kr and Xe, cal¬ 
culated using the full many-body theory (Dyson incident 
positron and full annihilation vertex) for s-wave incident 
positron. The corresponding tabulated values are given 
in Appendix B. For each of the noble gases shown, the 
valence np-subshell is the main component of the total 
7 -spectra at small Doppler energy shifts (< 2keV). This 
subshell is the most diffuse of the ground states. By con¬ 
siderations of the uncertainty principle, it therefore con¬ 
tributes to the spectrum mainly at low momenta (small 
energy shifts). Similarly, at small energy shifts, the mag¬ 
nitude of the valence shells is much greater than that of 
the subvalence or deeper lying shells, e.g., the valence 
shell of Kr is more than two orders of magnitude greater 
than that of the subvalence shell, which is a further three 
orders of magnitude greater than the next most inner 
shell. The need for the positron to overcome an increas¬ 
ing nuclear repulsion as it penetrates closer to the core 
causes this reduction in the probability of annihilation. 

The characteristic shapes of the individual subshell 
spectra are governed by the properties of the annihilated 
electrons. As the angular momentum of the annihilated 
electron increases, the spectra become flatter at small 
energy shifts. Indeed, for an s-wave positron it is easy 
to show from the first-order result (of the form of Eq. 13) 
and the small argument limit of the Bessel functions that 
the spectra behave as w( 0 ) — w(e) oc |e | 2 ^" +2 for energy 
shifts |e| < 1. Other features include the multiple ‘bell’ 
shapes, or ‘shoulders’, evident in the individual spectra. 
Looking at the subvalence 3s subshell of Ar in particu¬ 
lar, three distinct bell shapes are visible: the first in the 
region 0-2 keV, the second, with a much larger width, in 
the region 2-8 keV, and the third in the region > 8 keV. 
The number of distinct bell shapes in a given spectra is 
determined by the number of non-zero regions of energy 
density in the electron’s radial wave function, n—i. They 
can be understood by appealing to the uncertainty prin¬ 
ciple. For example, the 3s radial wave function of Ar, 
which has 2 nodes, has three regions of non-zero wave 
function density, wave function density at small radii 
manifests in the high-momentum (energy) regions of the 
spectra. The regions of non zero wave function density at 
larger r will contribute to the small momentum (energy 
regions) of the spectra. A further feature is that at large 
Doppler shifts the shape of the spectrum of the valence 
np-subshell, e.g., for Ar, the spectrum of the 3 p valence 
subshell follows that of the subvalence 2 p subshell; the 
requirement of orthogonality forces the 3 p and 2 p wave 
functions to be similar at small radii. Accordingly, their 
spectra are similar are large energy shifts (momenta). 

Table II shows the contribution of individual subshells 
to the total spectrum as a fraction of unity. Values 
are given for both the simplest calculation that uses a 
HF incident positron with the IPA vertex, and for the 
full MBT calculation (full vertex with Dyson incident 
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FIG. 7. Annihilation 7 -spectra for individual subshells of Ar, Kr and Xe, calculated with Dyson positron wave function and 
the full annihilation vertex. Only the valence n, sub-valence (n — 1) and (n — 2) subshells are shown for Kr and Xe as the more 
inner shells contribute negligibly to the total. 


TABLE II. Calculated fractional contribution to annihilation 
rate for noble gases (thermal incident positron). 


Subshell Argon Krypton Xenon 




HF-IPA 

MBT 

HF-IPA 

MBT 

HF-IPA 

MBT 

np 


0.8075 

0.8778 

0.8032 

0.8848 

0.7907 

0.8831 

ns 


0.1810 

0.1167 

0.1640 

0.998[-l] 

0.1616 

0.946[-l] 

(n - 

1 )d 



0.218[-1] 

0.107[-1] 

0.355[-l] 

0.175[-1] 

(n - 

1 )P 

0.830[-2] 

0.405[-2] 

0.852 [-2] 

0.365 [-2] 

0.973[-2] 

0.394[-2] 

(n - 

l)s 

0.314[-2] 

0.146[-2] 

0.244[-2] 

0.100 [-2] 

0.241 [-2] 0.926[-3] 

Total 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

(n - 

1) tot. 

0.011 

0.0055 

0.0328 

0.0154 

0.0476 

0.0224 

(n - 

1) tot. a 

<0.002 


0.013 


0.024 



a Calculated in Ref. [50] using semi-empirical fit. 


positron wave function). For all of the noble gases, the 
inclusion of the electron positron correlations result in an 
increased fraction of annihilation on the valence shells, 
compared to the core. In fact, the HF-IPA calculation 
overestimates the total core contribution ‘(n — 1 ) tot.’ by 
a factor of ~ 2. This is a result of the varying strength 
of the state dependent enhancement, which is stronger 
for the valence than for the inner shells. Also shown in 
the table are the results of the semi-empirical fitting of 
[50], in which the fractions were calculated by adjust¬ 
ing the ratio of HF-IPA inner to valence shell spectra in 
the total. For Kr and Xe this simplistic method gives 
fractional contributions that are of similar order to the 
parameterless MBT results, owing to the robustness of 
the shape of the spectra to the effect of correlations as 
discussed above (see Fig. 6 ). This method cannot how¬ 
ever, produce individual spectra of accurate magnitude 
like those of Fig. 7. Such spectra are required to calculate 
state-dependent 7 -spectra enhancement factors. 



FIG. 8 . Annihilation 7 -spectra of Argon. MBT calculated 
inner- (chain curve) and valence (dashed curve) shell contri¬ 
butions to the total (solid) are compared with experiment 
(circles) [50]. Also shown is the result of an IPA (zeroth- 
order vertex) HF calculation (dots). All calculated spectra 
have been scaled to experiment at zero energy shift. 


2. Comparison with experimental y-spectra: the 
core-electron contribution to the spectra 

Fig. 8 , 9 and 10 show the calculated spectra for Ar, Kr 
and Xe, after convolving with the experimental detector 
resolution function [50] 

/(e)=Wexp (iAk;) ■ (30) 
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FIG. 9. Annihilation 7 -spectra of Krypton (legend as de¬ 
scribed in caption of Fig. 8 ). 



FIG. 10. Annihilation 7 -spectra of Xenon (legend as de¬ 
scribed in caption of Fig. 8 ). 


where a = l/(41n2) 1 / 2 , TV is a normalization constant 
and AE’det = 1.16 keV is the width of the detector re¬ 
sponse, and scaling to the experimental spectra at zero 
energy shift. Specifically, the figures show the individ¬ 
ual contributions of the core and valence shells to the 
total spectra, calculated with the full MBT, the exper¬ 
imental results with associated errors, and the result of 
a simple HF-IPA calculation. In the experimental re¬ 


sults, the error bars represent the expected statistical 
variation in spectral amplitude due to counting statis¬ 
tics. Only the positive energy shifts of the symmetric 
7 -ray line is shown, since this region has greater statis¬ 
tical significance due to the reduced Compton scattering 
component [50]. In each figure the MBT calculated ‘va¬ 
lence’ spectra are the sum of the spectra of the ns and 
np subshells. As was seen in Fig. 7, the contribution of 
the (n — 2 ) and more inner shells are negligible com¬ 
pared to the subvalence shell. In Figs. 8 , 9 and 10 the 
‘inner’ spectra are therefore calculated as just the sum 
of the spectra for the individual (n — 1 ) subvalence sub¬ 
shells: (n — l)s + (n — l)p + (n — 1 )d, or just 2 s + 2 p for 
Ar. The relatively diffuse valence electrons, which have a 
narrow momentum distribution, dominate the spectrum 
at small energy shifts. This is clear from the figures, 
which show that the spectrum calculated from the va¬ 
lence shells alone agrees well with experiment at Doppler 
energy shifts < 2keV, but is insufficient at larger en¬ 
ergy shifts. In contrast, the core electrons are spatially 
compact, and therefore have a broader momentum dis¬ 
tribution compared with the diffuse valence electrons. 
They contribute markedly to the high-energy ‘wings’ of 
the spectrum, and their importance in bringing the total 
inline with experiment is evident for all three noble gases 
shown. For Ar, the inclusion of the core brings theory 
into excellent agreement with experiment over the full 
range of energy shifts. For Kr the valence shell con¬ 
tribution alone severely underestimates the total spec¬ 
trum, and the core contribution is of even greater impor¬ 
tance than in Ar, dominating over the valence at energies 
> 5 keV. Again the overall agreement of the spectra with 
experiment is excellent, although at higher energies the 
theoretical total slightly underestimates the experimen¬ 
tal spectrum. In Xe, the valence shell again severely un¬ 
derestimates the total spectrum. The core contribution 
dominates over the valence at energies > 4.5 keV, and is 
again seen to be vital in achieving agreement with exper¬ 
iment. Similar to the case of Kr, however, the theoret¬ 
ical spectrum slightly underestimates the experimental 
one at large energy shifts, and is ‘steeper’ than experi¬ 
ment. Since the theoretical results have been scaled to 
experiment, this suggests that the ratio of small to large 
momentum components is overestimated. 

One possible reason for the slight discrepancy is the ne¬ 
glect of relativistic effects on the electron wave functions. 
Indeed, Cook et al. [76] have shown that the momentum 
distributions of the valence Xe subshells show clear and 
direct relativistic effects manifestations of relativistic ef¬ 
fects. In general, the relativistic core orbitals will be 
more spatially compact than their non-relativistic coun¬ 
terparts and will have a correspondingly larger momen¬ 
tum density [77],[78]. Compared to the non-relativistic 
spectra, the high-energy regions of relativistic spectra 
should therefore be increased w.r.t. the low-energy re¬ 
gions, corresponding to an overall flattening. Such a flat¬ 
tening was evident in investigations of the Hartree-Fock 
and relativistic Dirac-Hartree-Fock (DHF) Compton pro- 
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files of the noble gases [79]. However, in that study the 
HF and DHF profiles were found to be similar at the 
momenta corresponding to the high-energy ‘wings’ of the 
Doppler spectrum and suggest that relativistic effects are 
not significant. In addition, since relativistic effects on 
the spectra scale as Z 2 [78], one would expect to see a 
larger discrepancy in the spectra of Xe (Z=54) than of 
Kr (Z=36). In contrast, we observe as similar discrep¬ 
ancy in both. A conclusive statement would require a 
full relativistic many-body calculation of the annihilation 
spectra. 

One further possible source of the discrepancy arises 
from the fact that, in order to be completely accurate, 
one should perform a Maxwellian average of the mag¬ 
nitude of the raw 7 -spectra over the incident positron 
momentum k, i.e., one should calculate the thermal av¬ 
eraged spectra given by 

{w{e)) k = (P/2n) 3/2 J w k (e)e~i 0k2 A-irk 2 dk, (31) 


IPA calculations to obtain accurate annihilation spectra. 

In general, one can define the momentum-dependent 
enhancement factor by approximating the non-local an¬ 
nihilation amplitude by the local, modified IPA form (c/. 
Eqns. (15) and (16)) 

Aie(P) ~ J e~ lF r tp e {r)ip n (r)y / j ne (r 1 P)dr, (32) 

where yj'fne is the state-dependent enhancement factor, 
usually interpreted as modifying the electron density at 
the positron (the square root is for later convenience). 
The calculation of 7 from many-body theory necessitates 
the reduction of our non-local annihilation amplitude 
to the IPA form of Eq. 32. Formally, using Eq. 16 and 
Eq. 32, the true enhancement factor is given by 

/ - t>n - 1 , If A P (r;rir 2 )V> e (ri)<p„(r 2 )dridr 2 
V7„«(r,F)-l+ fc(r)<ft,(r) 

(33) 


where j3 = 1 /ksT and ks is the Boltzman constant. 
This is likely to be broader than the spectra calculated 
at k = 0.04 a.u. 

Regardless, the discrepancy between the theoretical to¬ 
tal and the experiment at the high energy shifts in Kr and 
Xe is small, and despite it the MBT results show unequiv¬ 
ocally that: ( 1 ) the valence shells alone are insufficient to 
describe the total 7 -spectra of Ar, Kr and Xe, and that; 
( 2 ) there is a significant contribution from the relatively 
broad core spectrum to the total, and in addition, that 
this contribution is crucial to obtain agreement with the 
experimental spectra. 

B. Vertex enhancement factors for annihilation 
7-spectra 

It should by now be clear that the correlation correc¬ 
tions to the annihilation amplitude must be incorporated 
in an accurate calculation of the annihilation 7 -spectra 
of the valence and core subshells. The full many-body 
theory calculation is, however, non-trivial, and its ap¬ 
plication numerically intensive. To interpret the results 
of, e.g., PAES experiments, one needs some way of pre¬ 
dicting accurate core and valence annihilation spectra. 
Here, we appeal to the predilection of the literature and 
use the many-body theory to calculate enhancement fac¬ 
tors. These factors quantify the enhancement that results 
from the non-local corrections to the annihilation ver¬ 
tex. Such factors are common in positron-atom, positro- 
nium [11, 80], positron-charged-ion [45, 81] and positron- 
condensed-matter studies [2, 15, 26, 29, 30, 33, 34, 82, 83]. 
They are, however, usually calculated in some phe¬ 
nomenological way. They can be determined from an 
ab initio viewpoint from many-body theory, simply by 
calculating the required quantities in different approxi¬ 
mations to the vertex, and can be incorporated in simple 


The presence of nodes in the wave functions of the de¬ 
nominator, however, render this quantity of limited use 
and we must opt for a more pragmatic approach. 

The enhancement is in general dependent on the two -7 
momentum P [7]. However, we will see that to all inten¬ 
sive purposes this momentum dependence does not carry 
through to the annihilation 7 -spectra, at least for the core 
levels. It is therefore instructive to define a more sim¬ 
plistic state-dependent enhancement factor as the ratio 
of the calculated true partial annihilation rate (evaluated 
through Eq. 9) to that calculated using the zeroth-order 
(IPA) vetex; 

7 (o+i+r) 

_ _ /nA\ 


where the superscript denotes the vertex order (see 
Fig. 2) and n labels the electron subshell. This parame¬ 
terization is commonly used to predict true annihilation 
rates and 7 -spectra (see, e.g., [15, 26, 33, 34, 82, 83]). 
The true spectra for annihilation on a given subshell is 
then to be approximated by 

Wn(e) ~ 7nw[ 0) (e). (35) 

This procedure therefore assumes the vertex enhance¬ 
ment is independent of the Doppler-energy shift. 
Table III shows the values of 7 for the valence n and 
subvalence (n — 1) shells of Ar, Kr and Xe, calculated 
from the MBT using both the HF and Dyson incident 
positron. For the valence electrons, the vertex enhance¬ 
ment is dramatic, e.g., it produces a five-times increase 
in the spectral magnitude of the diffuse p-subshell of Ar, 
and a nine-times increase for that of Xe (for the HF 
positron). For the subvalence (core) subshells the factors 
satisfy 1< 7 CO re <7vaiencej as one should expect. They 
range from 1.35 (Ar 2s) to 1.87 (Xe 4d). These values 
are appreciably different from unity, and highlight the 
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TABLE III. Many-body theory calculated vertex enhance¬ 
ment factors 7 n(&Th) for noble gas subshells, calculated for 
s-wave thermal (fc = 0.04 a.u.) Hartree-Fock (HF) and Dyson 
(Dy) incident positron. 


Argon Krypton Xenon 


Subshell 

HF 

Dy 

HF 

Dy 

HF 

Dy 

np 

5.19 

4.42 

6.63 

5.23 

9.26 

6.57 

ns 

2.51 

2.39 

2.77 

2.57 

3.36 

2.99 

(n — 1 )d 



1.67 

1.62 

1.96 

1.87 

(n - l)p 

1.43 

1.42 

1.41 

1.38 

1.52 

1.48 

(n — l)s 

1.35 

1.34 

1.34 

1.32 

1.42 

1.39 

n tot. 

4.70 

4.02 

5.97 

4.73 

8.27 

5.89 

(n — 1 ) tot. 

1.41 

1.40 

1.56 

1.53 

1.85 

1.76 



FIG. 11. 7 -spectra for individual subshells of Krypton, calcu¬ 
lated using the full many-body theory (dashed and solid lines 
for core and valence subshells respectively) and reconstructed 
from the IPA spectra (for Dyson positron) using the enhance¬ 
ment factor 7 n in accordance with Eq. 35 (crosses for core and 
plus symbols for valence). 


necessary inclusion of correlations in the accurate calcu¬ 
lation of the core, as well as the valence, 7 -spectra. 

Reconstruction of the true 7 -spectra using these en¬ 
hancement factors via Eq. 35 is demonstrated in Fig. 11 
for the valence and subvalence subshells of Kr. Similar 
results are obtained for Ar and Xe [61]. For the sub¬ 
valence core (and the valence 4s) subshells, the energy- 
independent enhancement factors are seen to successfully 
reconstruct the true spectra. However, the reconstructed 
spectrum of the valence 4p-subshell overestimates the full 
MBT result at large Doppler shifts and underestimates 
at small ones. This is a result of the fact that both 



FIG. 12. Annihilation rate based enhancement factors 7ne 
for Ar (circles), Kr (squares) and Xe (diamonds), calculated 
using s-wave incident HF (solid symbols) and Dyson (unfilled 
symbols) positron. Dashed line is the fit (36) of 7 n i for atoms 
obtained using the static HF positron wave function (A = 42.0 
eV, B = 24.9 eV, /3 = 2.54), and the solid line is that for the 
Dyson positron wavefunction (A = 35.7 eV, B = 22.7 eV, /3= 
2.15). 


non-zero wave function density regions (‘lobes’) in the 
4 p radial wave function are scaled by the energy- and 
momentum- independent 7 factor equivalently. Since the 
core enhancement is smaller than the valence, we expect 
the lobe at smaller radial distances, which contributes to 
the large energy (momentum) regions of the spectrum, to 
scale by less than the lobe at larger r, which contributes 
to the small energy regions of the spectrum. 

In spite of these complications for the outermost 
valence-shells, the simple enhancement factors work re¬ 
markably well for the core shells. In Fig. 12 we have plot¬ 
ted the enhancement factors of Table III against the in¬ 
verse ionization energy of the subshell. The enhancement 
factor decreases in magnitude as the binding energy in¬ 
creases. 

As we have seen, the vertex enhancement to the core 
spectra is dominated by the first-order correction; higher 
order corrections result in a minimal amount of further 
enhancement (see Fig. 5 and Table I). One can estimate 
the strength of the short-range correlational enhance¬ 
ment from the first-order correction as the ratio Vint/A E, 
where V] n t is the electron-positron Coulomb interaction 
potential, and A E is a measure of the response of the 
electrons to the perturbation (of the order of the respec¬ 
tive electron’s binding energy). For positron annihila¬ 
tion on H-like ions, this strength scales as 1 /Z, where Z 
is the nuclear charge [7]. For annihilation on the core 
electrons of many-electron atoms the strength can be 
thought to scale analogously as 1 /VTe, where Ik is the 
ionization potential of the subshell. One also arrives at a 
1 /VTe dependence from simplifications to a calculation 
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FIG. 13. Enhancement factors (36) calculated using static 
HF (open symbols) and Dyson (solid symbols) positron states, 
for the s (triangles), p (squares) and d (circles) valence and 
core orbitals in Ar, Kr, and Xe; Is orbitals hydrogen (upside- 
down triangles) [7]; and hydrogenlike ions (plus signs) [7]. 
Dashed line is the fit (36) of 7 n i for atoms obtained using the 
static HF positron wave function (A = 42.0 eV, B = 24.9 eV, 
/3 = 2.54), and the solid line is that for the Dyson positron 
wavefunction (A = 35.7 eV, B = 22.7 eV, /3= 2.15). 

that explicitly includes two-particle Green’s function of 
the electron-positron pair in the annihilation vertex [61]. 
In light of this, we have fitted the function 



to the calculated enhancement factors of Table III, where 
A , B and /? are constants. The first term describes the 
enhancement of the core electrons and is physically mo¬ 
tivated as discussed above. The second term is phe¬ 
nomenological and is needed to describe the enhance¬ 
ment of the valence electrons, for which higher-order 
corrections are of as equal importance as the first-order 
one. The solid line in Fig. 12 is Eq. 36 fitted to the 
HF enhancement factors yielding the values A = 36 eV 
(1.32a.u.), B = 27.0eV (0.99a.u.), and /3 = 2.2. The 
dashed line is a fit of the function to the enhancement 
factors calculated with the Dyson incident positron wave 
function yielding the values A = 43.6eV (1.60a.u.), 
B = 20.5eV (0.75a.u.), and /3 = 2.6. 

Finally, Fig. 13 shows the values of the enhancement 
factors plotted against the ionization energy of the sub¬ 
shell, as well as the functional form Eq. 36 extrapolated 
to large Ie • The 1 /\/Te dependence for the core shells is 
evident. We also show the enhancements calculated for 
annihilation on Hydrogen and the H-like ions He + , Li 2+ , 
B 4+ and F 8+ . Given the inherent difference between the 
positive ions and the neutral systems, the ions follow the 
trend surprisingly well. For positrons annihilating in con¬ 
densed matter systems, the longer-range positron-atom 


correlations included in the Dyson positron quasiparti¬ 
cle wave function average to zero. In this case the HF 
vertex enhancement factors, which neglect these longer- 
range type correlations in the positron wave function, 
would be expected to be more appropriate. 


V. CONCLUSION AND OUTLOOK 

7 -spectra for positron annihilation on the core and va¬ 
lence subshells of the noble gases Ar, Kr and Xe has been 
calculated using a many-body theory framework. We 
have shown that the short-range electron-positron corre¬ 
lations, described through corrections to the zeroth-order 
(independent-particle-model) annihilation vertex, signifi¬ 
cantly enhance the 7 -spectra magnitude and annihilation 
rate parameter Z e ff. This enhancement is stronger for the 
valence shells, but it is still significant for the core shells 
(e.g., 7 ~ 1.9 for the 4 d Xe subshell). Despite the small 
probability of annihilation on the core subshells, we have 
shown that their contribution to the total annihilation 7 - 
ray spectra is vital for accurate description of experiment 
across the full range of Doppler-shift energies, and partic¬ 
ularly in the high-energy ‘wings’. We also use the atomic 
MBT to calculate true vertex enhancement factors that 
take account of the non-local corrections to the annihila¬ 
tion vertex. We have shown that for a given subshell n£, 
these factors can be described by a simple and physically 
motivated function of the subshell ionization energy I n E- 
Int - 1 = \JA/I nf , + ( B/I n tY , where A, B and /3 are 
positive constants. We have demonstrated the successful 
reconstruction of the core spectra using these enhance¬ 
ment factors to augment an independent-particle-model 
calculation. We suggest that this formula can be used 
to calculate the 7 -ray spectra for positron annihilation 
on the core electrons of atoms across the periodic table, 
and for the atomic-like core electrons of condensed mat¬ 
ter systems. This work, taken together with Refs. [5] and 
[ 10 ], gives a (near) complete description of the positron 
noble gas system. 

Looking to the future, there is still much work to be 
done to advance the theory and application of positron 
annihilation 7 -spectra. First, development and imple¬ 
mentation of a fully relativistic many-body theory would 
shed light on the discrepancies between theory and ex¬ 
periment in the high-energy ‘wings’ of the 7 -spectra of Kr 
and Xe. Moreover, accurate predictions of positron bind¬ 
ing and annihilation in complex molecules are required to 
develop new methods of vibrational spectroscopy [84] , to 
help explain the origin of the strong 7 -ray signal from 
the galactic centre [85], and are essential to realise spec¬ 
troscopic PET scanning [ 86 ] and advance antimatter- 
matter chemistry. Crucial in all cases is the ability of 
theory to accurately calculate the response of the atomic 
and molecular structure to the positron. For molecules, 
(crude) preliminary calculations of 7 -spectra for annihi¬ 
lation on small molecules, e.g., methane and its fluoro- 
substitutes, highlighted the strong effect of correlations 
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and the possibility of using positron annihilation to probe 
electron momentum densities [ 8 , 9]. With some develop¬ 
ments in the numerics, in particular an implementation 
of a Gaussian basis representation, accurate calculations 
for 7 -spectra on molecules could be calculated with the 
existing theoretical framework. 


In this expression the reduced delta-function element (c/. 
Eq. 13) is defined by 

Wa|| 2 ,i> = v^iM4] ^ ^ ^ 

x J j\(.Pr)Pe l t 1 i.r)P e2 t 2 {r)dr. 


Appendix A: Evaluation of y-spectra from angular 
decomposition of the annihilation amplitude 

The diagrammatic expansion of the annihilation ampli¬ 
tude has been given in Fig. 2. Its corresponding analytic 
form was given in Eq. 16. It is a sum of three terms, 
each of which is in one-to-one correspondance with the 
diagrams in Fig. 2. Writing the amplitude in angular 
decomposed form (making use of graphical angular mo¬ 
menta techniques [87]) and summing over the magnetic 
quantum numbers of the subshell m n , one can evaluate 
the y-spectrum as 

a p oo n+e 

*”' (e) = ai» E \Al2(P)\ 1 2 3 4 5 6 7 8 9 10 11 p<ip, (ad 

c A—|n—e| 

where 

A£(P) = (P\\®^\\ne), (A2) 

and the reduced vertex is given by 


where j\ is a spherical Bessel function, and we use the 
notation [f] = 2£ + 1. The reduced Coulomb matrix ele¬ 
ment for a particle pair propagating with coupled angular 
momenta A is 

(3,4||I/( a )|| 2,1) = ^(-1) a +'| * |* j4| {3)4 ||^|| 2)1)} 

(A5) 

where the reduced Coulomb matrix 

(3 ; 4||VJ||2,1) = 

/MMMM {ti £ *3 ) f 4 t h \ 

V 4tt ^oooyyoooy (A6) 

r r e 

x / Pe 3 e 3 (ri)Pe i eA r 2)-^lPe 2 i2{ r 2)Pe 1 eA r l) dr l dr 2- 
J r> 

Appendix B: Tabulated annihilation y-spectra for 
the noble gases Ar, Kr and Xe 



+ £ 


E ^ __ y- frJIM(^ll^ (A) 

E + e n — + e v 

V'” (A3) 

Aa | !a<- 2^2) (^2M2 I |r (A:> I Ia^I^i) (^lA^l 1 1 W (A) 

(E + £ n — £ lll + £ Vl )(E + £ n — £ M2 + e V2 ) 


ACKNOWLEDGMENTS 

We thank Prof. C. M. Surko for valuable discussions. 


[1] K. Iwata, R. G. Greaves, and C. M. Surko, Phys. Rev. 
A 55, 3586 (1997). 

[2] M. J. Puska and R. M. Nieminen, Rev. Mod. Phys. 66, 
841 (1994). 

[3] Core electrons, likewise, tunnel through the atomic bind¬ 
ing potential towards the positron, and their wave func¬ 
tions are affected by “exchange-assisted tunneling” [88]. 

[4] K. Iwata, R. G. Greaves, T. J. Murphy, M. D. Tinkle, 
and C. M. Surko, Phys. Rev. A 51, 473 (1995). 

[5] D. G. Green, J. A. Ludlow, and G. F. Gribakin, Phys. 
Rev. A 90, 032712 (2014). 

[6] L. J. M. Dunlop and G. F. Gribakin, J. Phys. B 39, 1647 
(2006). 

[7] D. G. Green and G. F. Gribakin, Phys. Rev. A 88, 032708 
(2013). 

[8] D. G. Green, S. Saha, F. Wang, G. F. Gribakin, and 
C. M. Surko, Mat. Sci. Forum 666, 21 (2010). 

[9] D. G. Green, S. Saha, F. Wang, G. F. Gribakin, and 
C. M. Surko, New J. Phys. 14, 035021 (2012). 

[10] D. G. Green and G. F. Gribakin, Phys. Rev. Lett. 114, 
093201 (2015). 

[11] G. F. Gribakin and J. Ludlow, Phys. Rev. A 70, 032720 
(2004). 


[12] K. G. Lynn, J. R. MacDonald, R. A. Boie, L. C. Feld¬ 
man, J. D. Gabbe, M. F. Robbins, E. Bonderup, and 
J. Golovchenko, Phys. Rev. Lett. 38, 241 (1977). 

[13] P. Asoka-Kumar, M. Alatalo, V. Ghosh, A. Kruseman, 
B. Nielsen, and K. Lynn, Phys. Rev. Lett. 77, 2097 
(1996). 

[14] M. P. Petkov, M. H. Weber, K. G. Lynn, R. S. Crandal 1, 
and V. J. Ghosh, Phys. Rev. Lett. 82, 3819 (1999). 

[15] M. Alatalo, H. Kauppinen, K. Saarinen, M. J. Puska, 

J. Makinen, P. Hautojarvi, and R. M. Nieminen, Phys. 
Rev. B 51, 4176 (1995). 

[16] A. Eshed, S. Goktepeli, A. Koymen, S. Kim, W. Chen, 
D. O’Kelly, P. Sterne, and A. Weiss, Phys. Rev. Lett. 
89, 075503 (2002). 

[17] A. Weiss, R. Mayer, M. Jibaly, C. Lei, D. Mehl, and 

K. G. Lynn, Phys. Rev. Lett. 61, 2245 (1988). 

[18] T. Ohdaira, R. Suzuki, T. Mikado, H. Ohgaki, M. Chi- 
waki, and T. Yamazaki, Appl. Surf. Sci. 116, 177 (1997). 

[19] J. Mayer, K. Schreckenbach, and C. Hugenschmidt, 
Phys. Status Solidi (C) 4, 3928 (2007). 

[20] A. H. Weiss, N. G. Fazleev, M. P. Nadesalingam, 
S. Mukherjee, S. Xie, J. Zhu, and B. R. Davis, Radiat. 
Phys. Chem. 76, 285 (2007). 









18 


TABLE IV. Annihilation 7 -spectra w(e) for valence and subvalence shells of the noble gases Ar, Kr and Xe. Results quoted 
are that of the full MBT calculation, i.e., fully-correlated Dyson incident positron and full annihilation vertex. Square brackets 
denote powers of ten. Spectra in different approximations to the incident positron wave function (HF or Dyson) and different 
orders of the annihilation vertex can be requested from the authors. 


e (keV) 

2 s 

Argon 

2 p 3s 

3 p 

Krypton 

3s 3p 3 d 

4s 

4 p 

Xenon 

4s 4 p 4 d 

5s 

5 p 

0.00 

7.28[-3] 

1.17 

1.59 

8.96 

1.60[-2] 3.57[-2] 8 .93[-2] 

4.05 

2.77[1] 

9.15[-2] 2.44J-1] 9.38[-l] 

2.14[1] 

1.53[2] 

0.25 

7.24 [-3] 

1.17 

1.51 

8.92 

1.59[-2] 3.57[-2] 8 .93[-2] 

3.79 

2.74[1] 

9.02[-2] 2.44J-1] 9.37[-l] 

1.98[1] 

1.51 [2] 

0.50 

7.10[-3] 

1.17 

1.31 

8.62 

1.54[-2] 3.56[-2] 8 .93[-2] 

3.17 

2.60[1] 

8.63[-2] 2.44[-l] 9.37[-l] 

1.56[1] 

1.41 [2] 

0.75 

6.88 [-3] 

1.17 

1.03 

7.83 

1.48[-2] 3.56[-2] 8.92[-2] 

2.38 

2.27[1] 

8.01[-2] 2.43[-l] 9.37[-l] 

1.06[1] 

1.17[2] 

1.00 

6.59[-3] 

1.16 

7.33[-l] 

6.51 

1.39[-2] 3.54[-2] 8.92[-2] 

1.61 

1.78[1] 

7.22[-2] 2.40[-l] 9.35[-l] 

6.11 

8.29[1] 

1.25 

6.23[-3] 

1.16 

4.78[-l] 

4.94 

1.28[-2] 3.51 [-2] 8.90[-2] 9.75[-l] 

1.26[1] 

6.32[-2] 2.35[-l] 9.30[-l] 

3.08 

5.02[1] 

1.50 

5.82[-3] 

1.15 

2.85[-l] 

3.44 

1.16[-2] 3.46[-2] 8.87[-2] 5.32[-l] 

7.93 

5.36[-2] 2.26[-l] 9.18[-1] 

1.45 

2.61 [ 1 ] 

1.75 

5.36[-3] 

1.14 

1.57[-1] 

2.21 

1.04[-2] 3.38[-2] 8.80[-2] 2.70[-l] 

4.49 

4.41 [-2] 2.14J-1] 8.97[-l] 7.96[-l] 

1.18[1] 

2.00 

4.89[-3] 

1.12 

8.29 [-2] 

1.33 

9.07[-3] 3.27[-2] 8.69[-2] 1.42[-1] 

2.29 

3.52[-2] 1.99[-1] 8.63[-l] 6.24[-l] 

4.69 

2.25 

4.40 [-3] 

1.09 

4.53[-2] 7.57[-l] 

7.79[-3] 3.14[-2] 8.51[-2] 9.47[-2] 

1.05 

2.73[-2] 1.81 [-1] 8.17[-1] 6.02[-l] 
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